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Quantum walks have emerged as an interesting alternative to the usual circuit model for quantum 
computing. While still universal for quantum computing, the quantum walk model has very different 
physical requirements, which lends itself more naturally to some physical implementations, such as 
linear optics. Numerous authors have considered walks with one or two walkers, on one dimensional 
graphs, and several experimental demonstrations have been performed. In this paper we discuss 
generalizing the model of discrete time quantum walks to the case of an arbitrary number of walkers 
acting on arbitrary graph structures. We present a formalism which allows for analysis of such 
situations, and several example scenarios for how our techniques can be applied. We consider the 
most important features of quantum walks - measurement, distinguishability, characterization, and 
the distinction between classical and quantum interference. We also discuss the potential for physical 
implementation in the context of linear optics, which is of relevance to present day experiments. 



I. INTRODUCTION 

Quantum walks [1 offer an alternate approach to im- 
plementing quantum algorithms [2] compared to the typi- 
cal circuit [3^ or cluster state [11 15] models for quantum al- 
gorithms. Indeed, quantum walks have been successfully 
applied to database search [6] and graph isomorphism 
|7 problems. More recently, the universality of quantum 
walks has been shown in the case of continuous [8 and 
discrete [9 systems. A review of both continuous and 
discrete walks can be found in Ref. [10 . 

For more than a decade the quantum walk was a very 
fruitful yet theoretical concept. However, experiments 
with single optically trapped atoms [11 have successfuly 
demonstrated that quantum walks can be realized in the 
laboratory. A quantum walk on trapped ions ^\ has 
been realized. Various successful realizations of a quan- 
tum walk based on quantum optics implementation fol- 
lowed. 

Indeed, photonic discrete time quantum walks were ex- 
perimentally demonstrated by Schreiber et al. [13] us- 
ing weak coherent light and by Broome et. al [T4^ using 
single photons produced via parametric down-conversion 
(PDC). The first experiments investigating the impact 
of correlations in quantum walks with two walkers were 
done in the continuous time model: Bromberg et. al 
[T5] showed correlations by classical intensity correlation 
measurements and Peruzzo et. al fW^ used coincidence 
measurements of two photons of a PDC source, both with 
integrated wave guide devices. 

Most of the studies to date have focused on walks 
with a single quantum particle. Nevertheless, some works 
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have considered the effect of multiple walkers on a line. 
Indeed, quantum walks involving more particles unlock 
additional possibilities which are unavailable in classical 
random walks, as the quantum particles can be entangled 
or they can be indistinguishable fermions or bosons. In 
Ref. [17 the situation with two perfectly correlated walk- 
ers is considered, while in Refs. [181 El the two walker 
situation is considered with separable and initially en- 
tangled walkers. The meeting problem in this model has 
been analyzed in Ref. [20 . Additionally, Ref. [21] con- 
siders many-body entanglement using one dimensional 
walks. 

In this paper we present the following: (1) a framework 
for multi-walker walks on arbitrary graph structures; (2) 
the difference between distinguishable and indistinguish- 
able walkers, and the issue of walker correlations; (3) the 
characterization of quantum walks; (4) the effects of mea- 
surement on multi- walker systems in the context of both 
distinguishable and indistinguishable walkers; and, (5) we 
turn our attention to the example of a two photon lin- 
ear walk with both distinguishable and indistinguishable 
walkers, which is of relevance to present day experiments. 
Our work builds on previous work on quantum walks by 
introducing techniques for modeling multiple walkers and 
arbitrary graph structures. Until now there has only been 
very limited consideration of multi-walker systems. We 
do this by introducing a walker operator technique for 
describing such systems. In the context of the multiple 
walker scenario we discuss the main important issues sur- 
rounding measurement and characterization, and present 
a detailed discussion of the special case of photonic im- 
plementation. 

The archetypal discrete time quantum walk consists of 
two building blocks - a coin operator and a step oper- 
ator. This terminology arises historically in accordance 
to models of discrete classical random walks. For consis- 
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tency we kept the formalism, although the two operators 
can be merged into a single evolution operator. We de- 
fine a basis of states of the form |x,c), where x denotes 
a position state of a walker^ and c denotes a value of the 
coin. The position state could literally refer to the phys- 
ical position of a particle, as in the case of most pho- 
tonic implementations, or it could refer more abstractly 
to the occupation of a level in some higher dimensional 
degree of freedom. The coin is essentially an ancillary 
parameter which is used by the step operator to decide 
how to propagate the walker. The simplest example of a 
quantum walk is the one dimensional case with a single 
walker, where the positions x represent positions along a 
one-dimensional lattice graph. A common choice of coin 
and step operators are, 

H\x,±l) = {\x,-l)±\x,l))/V2 

S\x,c) = |x + c,c) (1) 

where c G { — 1, 1} and x G {Z}. H is the Hadamard coin, 
and S is the nearest neighbor step operation. Thus, the 
coin operator affects only the coin state and leaves the 
position state unchanged, whereas the step operator de- 
termines the evolution of the position state, as a function 
of the value of the coin, leaving the coin state unchanged. 
This particular choice of coin and step operators can be 
optically implemented using 50/50 beamsplitters. After 
n rounds of the quantum walk, the evolution is of the 
form {SH)'^. At each step the reach of the position state 
expands outwards across the line. Additionally, because 
the position states propagate in both directions, inter- 
ference between the position states takes place. This is 
the defining feature that differentiates the quantum walk 
from a classical random walk [27] , 

In the above description of a quantum walk we have 
made a specific choice of coin and step operators. Of 
course, we have great fiexibility in choosing these oper- 
ators. In fact, the algorithm implemented by a quantum 
walk is characterized entirely by this choice. For example, 
we could use a biased coin, which can optically be imple- 
mented using a beamsplitter with refiectivity rj^ ^ 1/2, 
or a step operator with, for example, periodic bound- 
ary conditions, S\x^c) = | (x + c) mod A^, c) , where N is 
the number of distinct position states. Not only can we 
choose arbitrary operators, but, in general, the choice can 
vary at each step of the walk. In the following sections we 
will describe techniques for analyzing multi-walker walks 
on arbitrary graphs, and some of their properties. 



II. MULTI- WALKER QUANTUM WALKS ON 
ARBITRARY GRAPHS 

We now turn our attention to generalizing the quan- 
tum walk to the case of multiple walkers and arbitrary 
graph structures. This is motivated by the increasing 
complexity of the system as the number of walkers in- 
creases. Additionally, algorithms are emerging relying on 



multiple walkers [22] and higher dimensional graph struc- 
tures [23 . First let us develop a notation for representing 
multi-walker states. We begin by introducing the walker 
operator, w{x,cy, which creates a walker at position x 
with coin value c. We also define the empty walker state, 
|0), which is the state of the system without any walkers 
present. In an optical context the walker operator cor- 
responds directly to the photon creation operator, and 
the empty walker state is the vacuum state. The walker 
operator has the following effect, 

w{x,cy\0) = \x,c) (2) 

and all walker operators commute. An arbitrary initial 
state of a quantum walk system with multiple walkers 
can be represented as combinations of different walker 
operators. 

Next we introduce the walker graph, which captures the 
available position states, and the relationship between 
them. In general these graphs will be bidirectional, with 
weighted edges. Some examples are illustrated in Figs, [l] 
and[2j We define unitary operators A^^\ which act on the 
neighborhood of x. The neighborhood of position state i is 
denoted and consists of i and its immediate neighbors 
with outgoing directed edges. 

Now we define the coin and step operators as follows, 

S: w{xjy ^ w{j,xy (3) 

where A is unitary for all x. Note that the definition of 
the step operator in this case is not quite analogous to 
the simple example presented in Eq. [l] There the step op- 
erator updates the position while leaving the coin state 
unchanged, whereas in the above formalism the step op- 
erator is a permutation operator that changes both the 
coin and the position. This slight change in the formal- 
ism is necessary to ensure unitarity of the step operator. 
For example, if in analogy to Eq. [l] our step operator 
implemented the transformation w{x,j)^ w{j,j)\ the 
operation would not be reversible. 

The interpretation of the coin state varies at different 
points in the evolution of the walker. Prior to the coin 
operator, the coin value tells us the source of the walker 
(i.e. its previous position state), whereas after the coin 
operator it tells us the destination of the walker (i.e. its 
next position state). This is consistent with the defini- 
tion of the linear walk in Eq. [T] Thus the coin operator 
puts the walker into a superposition of destination states 
within the neighborhood, and the step operator updates 
the walker to the new destination while replacing the 
source state with the previous position state. In general, 
time-dependent coins are allowed, and the evolution of 
the system will be of the form Y\i^^ii where Ci is the 
coin operator at time step i. 

This formalism maps a position to a destination, 
with memory of the previous state. In general, arbitrary 
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graphs will have varying degree, i.e. different neighbor- 
hoods contain different numbers of vertices. Thus in gen- 
eral we cannot use a single coin to describe the dynamics 
of the entire graph. Instead, unique operators must 
be defined for every neighborhood. This differs from our 
initial example of a walk on a line, where the degree is 
constant and a single coin can characterize the entire sys- 
tem. 

Importantly, in the above formalism we have implicitly 
assumed that the coin is independent of the number of 
walkers. For example, in the optical context one might 
imagine that the coin is based on a cross-Kerr effect, 
in which case the coin is dependent on the number of 
walkers. We implicitly rule out these kinds of effects such 
that, for examp le, i n the optical case (to be discussed 
in detail in Sec. Ill), propagation is defined entirely by 



linear elements - beamsplitters and phase shifters. 

Using this formalism the walk is characterized com- 
pletely by (1) the initial state of the system, and (2) the 
unitary coin matrices of the graph. Additionally, the evo- 
lution of a particular walker is characterized entirely by 
its neighborhood in the graph. The evolution of the sys- 
tem, and the algorithm implemented by it, are specified 
solely by these parameters. 



A. Examples of graphs 

1. One- dimensional lattice graph 

As a first simple example we will step backwards to the 
simple example of walkers on an infinite one-dimensional 
lattice graph. This graph is already well understood, so 
we will simply demonstrate in this section that the gen- 
eral formalism we have presented reduces to the dynamics 
presented earlier in Eq. [l] However our formalism accom- 
modates for arbitrary number of walkers rather than just 
one. The position graph for the linear walk is shown in 
Fig.[l] 

FIG. 1: Position graph for a one dimensional linear walk. 

To model this graph in our general formalism, we only 
have to chose appropriate operators A^^^ . For this graph 
we choose the matrix elements of the operators as follows, 

nx{l),nx{l) 



Ax) 

V(l),n.(2) 



A 



ix) 

na;(2),n,c(l) 



= 1/V2 Vx 
= 1/^2 Vx 
= I/V2 \/ X 



Ax) 

^n,(2),n^(2) 



1/V2 M X 



(4) 



operator, acting on the subspace spanned by n^^ in anal- 
ogy to the linear system presented earlier. Note that we 
have made the coin operator the same for all values of 
X, which is legitimate when the line is infinitely long (i.e. 
without boundary conditions). 



2. Two dimensional lattice graph 

Next we turn our attention to a simple two dimensional 
example - a two dimensional regular lattice. The position 
graph is shown in Fig. [2] In this example we will let the 
operators A(^) be the same for all values of x. Again, in 
general this would not apply with boundary conditions, 
but on an infinite lattice this is a valid assumption. 




FIG. 2: Position graph for a walk on a two dimensional lattice. 

An example of a coin on this graph is the four level 
Or over coin, given by 
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where nx{i) is the zth neighbor of x - in this case there 
are only two. Thus the operator A^^^ is the Hadamard 



which can be generalized to higher dimensions. This coin 
is essentially a generalization of the Hadamard coin to 
a four level system [28 , and for a given position, defines 
the mapping between the old coin value and the new coin 
value. 



B. Distinguishable versus indistinguishable walkers 

In the above discussion we have implicitly assumed 
that walkers with the same position and coin values are 
indistinguishable, i.e. if we have two walkers in a given 
position and coin state, measurement cannot determine 
which walker is which. Should we want to model dis- 
tinguishable walkers we can introduce a new parameter 
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which identifies a given walker from the others. In this 
case we redefine our walker operator as 



(6) 



where n is a new parameter which identifies which walker 
the operator corresponds to. In an optical context this 
could, for example, identify in which spatio-temporal 
mode or polarization the photon is in. We now have the 
property that {0\w{x^ c^n)w{x^ c^my\0) = Sm,n- In the 
following sections we will primarily consider the situation 
of indistinguishable walkers, since this is the interesting 
case where the walkers interfere with one another. How- 



ever, in Sec. IIII CI we will consider the two cases in the 



optical context and examine the difference in the inter- 
ference patterns. 



C. Characterizing quantum walks 

A completed quantum walk, up to phases, can be fully 
characterized by a joint probability distribution (JPD) 
[18l[T9]. The JPD has dimension n, where n is the num- 
ber of walkers, each degree of freedom represents the po- 
sition of a given walker, and each element in the JPD 
is the joint (i.e. coincidence) probability of measuring a 
walker at each of the corresponding positions. We define 
Pii,i2,...,iri to be the probability that a walker is found at 
each of the positions n, ^2, • • • , (an n-fold coincidence 
probability) . Optically this would simply be implemented 
by placing photo-detectors at each position and measur- 
ing joint detection events. 

The number of terms in the JPD is 0{N^)^ where TV 
is the number of distinct position states. This implies 
that fully characterizing the statistics of a multi-walker 
quantum walk is infeasible with large numbers of walk- 
ers, necessitating the development of algorithm specific 
metrics for characterizing the system which do not re- 
quire exponentially large numbers of measurements. This 
scaling raises the interesting side-thought as to how the 
computational power of a quantum walk scales with the 
number of walkers - with a single walker the number of 
parameters characterizing the system is 0{N)^ whereas 
for n > 1 the number of parameters grows exponentially 
with n. The JPD can be used to derive any output statis- 
tics of interest. 

As an example, in the meeting problem we are inter- 
ested in the probability that all walkers are measured at 
the same position, which is given by tr(P). The meeting 
problem is one of the standard parameters discussed in 
connection with classical walks involving multiple par- 
ticles [24^. The meeting probability quantifies the pos- 
sibility the two walkers meet. In the quantum domain 
the meeting probability quantifies the analogous property 
but takes also into account additional quantum proper- 
ties of the walkers. Hence it is natural to expect quite a 
different behavior of the meeting probability when com- 
pared to the classical case. The effects of walk topology 



and indistinguishability of the walker can be nicely repre- 
sented. On the short-time scale, quantum features result 
in strong increase of the meeting probability for indistin- 
guishable bosons, an effect familiar to Hong-Ou-Mandel 
interference. Moreover, the differences in the meeting 
probability can be found also in the asymptotic regime. 
However, the quantum effects on the long-time scale are 
much more subtle, as discussed in detail in Ref. ^20^ . 



D. Measurement 

We now investigate the effects of measurement on a 
multi-walker system. We will focus on the case of two 
walkers, although this could be easily generalized to arbi- 
trary numbers of walkers. Consider a two- walker system 
on an arbitrary graph with the first walker initialized into 
position X and the second into position y. Then following 
n applications of coin and step operators, U = {SC)^, 
we obtain 

w{x)^w{y)^ ^ [Yl U.Mi)^) ( E UyMJ)^) (7) 

i 3 

where we have ignored coin values for simplicity. Note 
that although this state is written in a factorizeable form 
it is nonetheless entangled, since it cannot be rewritten 
in a form which is a product of the walkers at different 
positions. Next, suppose we perform a measurement and 
detect a single walker at position m. Then the state will 
be projected into. 



(8) 



Thus, upon measurement the factorizeable nature of 
the state is lost, and our resulting state contains cross- 
combinations from the different terms in the initially fac- 
torizeable state. Therefore the behavior of a mult i- walker 
system cannot, in general, be characterized in terms of 
single walker systems post-measurement. 

In an optical context this corresponds to preparing 
photons at positions x and allowing them to propa- 
gate through a linear optics network (beamsplitters and 
phase-shifters) characterized by the unitary matrix /7, 
and then measuring a single photon at position m. 

We could also characterize the system in terms of the 
coincidence probability, which, for two walkers, is given 
by [15^ 



y,n 



(9) 



Again, in an optical context this equates to measuring a 
photon at both positions m and n - a standard coinci- 
dence measurement. 

The above discussion applies for indistinguishable 
walkers. In the case of distinguishable walkers the state, 
upon measurement, is projected onto 



w{m, iyUx,r 



^2)^) 



(10) 
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where 1 and 2 denote the classicahy distinguishable walk- 
ers and we have measured walker 1. Similarly, the coin- 
cidence probability for distinguishable walkers is given 

by 



I Ux.mUy 



(11) 



where we have measured the first walker at position m 
and the second at n. This expression contains no interfer- 
ence terms, as expected. It is evident from Eq.[TT]that for 
distinguishable walkers the coincidence probability has 
an uncorrelated form, i.e. it can be represented as the 
product of two independent functions. Hence the coinci- 
dence statistics for distinguishable walkers can be char- 
acterized in terms of single walker walks. This will be 



discussed further in the optical context in Sec. |III C 



An interesting observation is that when both walkers 
are initialized into the same position state, x = the 
coincidence probability Pm,n is of the same form between 
the distinguishable and indistinguishable scenarios [29 . 
This implies that we cannot expect a quantum speedup 
over the single walker case by preparing multiple walkers 
into the same initial state. 



III. APPLICATION TO QUANTUM OPTICS: 
TWO PHOTON WALK ON A LINE 

A. The two photon beamsplitter coin 



sub-spaces: the zero photon block implements the trivial 
identity operation; the single photon block implements 
the usual Hadamard coin; and, the two photon block im- 
plements an unbalanced generalization of the Hadamard 
coin, with its top-left element zero, corresponding to 
Hong-Ou-Mandel (HOM) [25 interference between the 
incident photons. In the single walker case, the beam- 
splitter coin will consist of just the first two blocks. Im- 
portantly, the Hq and Hi blocks can be simulated using 
coherent light, whereas the H2 block cannot be, as it 
contains HOM interference terms. This differentiates the 
single and two walker scenarios. 

Experimentally, the beamsplitter simultaneously im- 
plements the coin and step operators. That is, it simul- 
taneously flips a quantum mechanical coin and subse- 
quently propagates the physical state spatially. A simple 
optical quantum walk using beamsplitter coins is illus- 
trated in Fig. [3) This configuration is ideally suited to 
implementation using integrated waveguides, as has been 
demonstrated in Ref. [H]. 



step 




We now turn our attention to the illustrative example 
of the two photon, one dimensional quantum walk. Sup- 
pose we have a system containing exactly two photons 
(walkers), and the propagation of the quantum walk is 
performed by beamsplitters. Then each input to a beam- 
splitter will have 0, 1 or 2 photons. We can define the 
basis states for the input and output of the beamsplitter 
as |00), |01), |10), |11), |20) and |02), where \mn) is the 
two mode state with m photons in the first mode and n 
photons in the second. For now we assume the incident 
photons are indistinguishable. Acting on this basis, the 
beamsplitter transformation can be represented in matrix 
form. 



Bo 




(12) 



This follows directly from the beamsplitter equations of 
motion, {a^ + b^)/V2 and {a^ - h^)/V2 

(using the asymmetric beamsplitter convention), where 
a and h denote the two input/output modes. The ma- 
trix is block diagonal, B2 = Hq i^i i^2 7 with each 
of the blocks representing the different photon number 



FIG. 3: A simple optical 3 level quantum walk on a line. 



Photon statistics for the two photon quantum 
walk 



We ran a simulation with both 4 levels and 7 levels 
in the tree, making 8 and 14 output modes respectively. 
In Figs. [4] and [5] we plot the probability amplitudes for 
the different terms in the output expression as a JPD. 
Let the numbers on the axes be denoted by m and n. 
Then the element {m, n} corresponds to the term with 
a photon at mode m and a photon at mode n. The ma- 
trix is symmetric since there is no way of distinguishing 
{m, n} from {n,m}. The diagonal terms correspond to 
the terms where both photons are in the same output 
mode, while the off-diagonal terms represent the terms 
where the two photons are in different output modes. 
This leads to a normalization condition for the JPD ^O] . 
This graph characterizes all the coincidence probabilities 
associated with the different terms in the output state. 
Note that all of the off-diagonal elements of the JPD can 
be experimentally obtained using coincidence measure- 
ments, which does not require number resolving detec- 
tion. However, the diagonal terms, in which both pho- 
tons appear at the same output (the meeting problem), 
inherently requires number resolving detection. 
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FIG. 4: (Color online). Two photon quantum walk with 4 levels, (left) Joint probability distribution with two indistinguishable 
walkers, (right) Joint probability distribution with two distinguishable walkers. Note that in the case of distinguishable walkers 
the coincidence matrix is uncorrelated and can be obtained by multiplying the single photon statistics along one axis by 
the single photon statistics along the other axis. In the case of indistinguishable walkers the matrix is correlated, which is 
characteristic of a quantum behavior. 



(a) (b) (c) 




First phoion First photon 



FIG. 5: (Color online). Two photon quantum walk with 7 levels, (a) Joint probability distribution with two indistinguishable 
walkers, (b) Joint probability distribution with two distinguishable walkers, (c) Single photon statistics for both distinguishable 
and indistinguishable walkers. Visually there is no discernable difference in (c) between the distinguishable probabilities (Pj;^'') 
and indistinguishable probabilities (P/^'*), although there is a small numerical difference Li = | \Pj;'^'' — Pi'^''\ = 0.04. As 
with the 4 level system, in the case of distinguishable walkers the plot is uncorrelated, whereas in the case on indistinguishable 
walkers it is correlated. 



Next we consider the probability of detecting one or 
more photons at the mth mode, i.e we are now consider- 
ing single click events rather than coincidence events. If 
we take the mth column from the previous matrix and 
sum its elements, we obtain the net probability of detect- 
ing a photon in mode m. This is shown in Fig.jsjc). This 
is analogous to the single walker case where we plot the 
probability distribution of the single photon statistics. 
Interestingly, the projection onto single click events in 
the case of indistinguishable and distinguishable walkers 
produce very similar statistics. The full matrix, however, 
contains much more information as it characterizes both 
single photon and coincidence statistics. 

The two photon JPD and its projection onto one di- 
mension can both be directly experimentally measured. 
The former can be measured by performing coincidence 
measurements, which tells us the probability associated 
with the respective terms, {m, n}. The projection can be 
measured using single (non-coincidence) measurements 
at the outputs. 



Notice that in the JPD of indistinguishable particles 
exist clear oscillations between different terms, which are 
distinct from the distinguishable case. This structure is 
lost upon projection into one dimension. In the case of 
the 4 level system, it is evident that the two photon in- 
terference is characterized by the suppression of the two 
diagonal peaks. In the case of the distinguishable photons 
the JPD is uncorrelated, whereas in the indistinguishable 
case it is not. 



An interesting question is, how do we define the 'out- 
put' of the walk? Is it given by the single photon or the 
coincidence statistics? If we define the output as being 
the single photon statistics then the single and two walker 
situations are very similar - they produce nearly identi- 
cal spectra. However, if we define the output as being the 
set of coincidence statistics, then we have a much richer 
data set, which does not have a counterpart in the single 
walker scenario. 
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C. The transition from distinguishable to 
indistinguishable walkers 



The analysis has, up until this point, considered the 
situation of two indistinguishable walkers. However, in 
an arbitrary quantum walk, distinguishable walkers are 
also possible with an appropriate choice of coin and step 
operators. This arises naturally with distinguishable pho- 
tons. We now turn our attention to comparing the situ- 
ations where the walkers are indistinguishable and dis- 
tinguishable, and examine the differences in the coinci- 
dence statistics in both scenarios. Let the input modes 
to the first beamsplitter be denoted and 1. Then, when 
the incident photons are in the same spatio-temporal 
mode, the input state can be expressed ajaj. On the 
other hand, when the incident photons are distinguish- 
able, we can express the input state as aJ^J, where b 
denotes a mode orthogonal to a. In the former case the 
incident photons will interfere within the quantum walk, 
whereas in the later case they will not and instead evolve 
completely independently of one another and no inter- 
ference will be observed. In the indistinguishable case we 
obtain JPDs given in Figs. [4]^left) andjsja), whereas in 
the distinguishable case the JPDs are illustrated in Figs. 
[Ij^right) and |5|b) [31]. Notice the differing interference 
patterns between the indistinguishable and distinguish- 
able walker scenarios [15 . However, upon projection onto 
single click events, both scenarios produce similar statis- 
tics. In the case of distinguishable walkers, the JPD is 
given by the product of the single click statistics along 
each axis, whereas in the case of indistinguishable walkers 
the JPD is correlated. 

In the case of distinguishable walkers, {m, n} and 
{n,m} are now distinct. This leads to a different nor- 
malization condition than for the indistinguishable case 

m- 

The two situations above, of indistinguishable and dis- 
tinguishable walkers, represent the two extreme cases. In 
the in-between scenario whith an arbitrary degree of in- 
terference, we can decompose the two photon state into 
overlapping and non-overlapping components [26] and 
represent the input state as aajaj + Vl — a'^alb\, where 
a G 3? quantifies to what extent the walkers are indistin- 
guishable or distinguishable. In this case, the probabili- 
ties in the JPD are simply given by a linear combination 
of the two extreme cases [33^ . 

Note that there are two ways of interpreting the pa- 
rameter a. On one hand, with complete control over the 
system, it can be used to 'tune' the extent to which the 
walkers are indistinguishable or distinguishable. On the 
other hand, it could be forced upon the experimenter as 
a result of mode-mismatch, time-jitter, or other forms of 
distinguishability between the incident walker photons. 



D. Photon correlation 

As noted above, the defining feature of walks with in- 
distinguishable walkers is the correlation that exist be- 
tween the photons in the system, which is not the case 
for distinguishable walkers. This arises as a result of in- 
terferometric suppression or enhancement of particular 
terms in the JPD. This was experimentally observed by 
Bromberg et. al in a continuous time quantum walk [15]. 
More formally we can define a correlation measure as fol- 
lows: (1) write the JPD in a matrix form, (2) diagonalize 
the JPD matrix with diagonal elements given by p^, (4) 
renormalize the elements such that 

Pi ^ ^ (13) 

gives a probability distribution, and (5) calculate the 
Shannon entropy of the resulting probability distribution 

H{P) = -J2Pi^ogPi. (14) 

i 

As an example we consider the 4 level walk presented 
earlier. For distinguishable walkers the JPD is uncorre- 
lated, so it is diagonalizable with a single non-zero coef- 
ficient, hence pi = 1 and H{P) = 0. On the other hand, 
for indistinguishable walkers H{P) is necessarily greater 
than zero. In this example we have H{P) = 1.74, a char- 
acteristic of photon correlations. 

E. Comparison with coherent state 
implementations 

We conclude our discussion of optical implementa- 
tions of quantum walks by comparing to implementa- 
tions where coherent states, rather than single photons, 
are employed. Such a system was recently demonstrated 
by Schreiber et. al p^. 

We show that an arbitrary quantum walk system, in 
which all walkers are initialized into the same position 
state, can be simulated with coherent light (i.e. both will 
yield identical output statistics). This is consistent with 
previous observations in the continuous time regime [15j. 
On the other hand, systems where walkers are initial- 
ized into different position states cannot, in general, be 
simulated by coherent light. 

The evolution of a system containing n steps, with p 
walkers all initialized into position x, is given by 

(f^4)^|o) = (E^fc4)'lo), (15) 

k 

where U = (SC)^. Let us now replace the p-photon Fock 

state with a coherent state \a) = e~^e^^i|0), also ini- 
tialized at position x. The evolution is given by 

= e-^e^/^^^^^^HO) (16) 
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The second line follows from 15 Performing a Taylor se- 
ries expansion we obtain 



e 2 e 



d ' k 



(17) 

If we now condition on p walkers {d = p) (in the optical 
context using photon number resolving photo-detection) 
we obtain the same state as in the Fock state example, 
up to a constant factor c 



.ML 



f(E«4)'io) 



(18) 



Note that the state is of a separable form. In our previ- 
ous discussions we observed that walks consisting of two 
walkers at different inputs generate non-separable coinci- 
dence statistics, which is impossible in this case. Thus, in 
general, coherent states can only simulate multiple walk- 
ers in special cases: (1) with distinguishable photons, 
which generate separable statistics, and (2) the exam- 
ple presented above where multiple walkers begin in the 
same initial state. Importantly, even with quantum mea- 
surements, such a system cannot mimic the statistics of 
a system with quantum states. 



IV. CONCLUSION 



Since c is constant it does not change the distribution of 
the single click or coincidence statistics and the distribu- 
tion reflects that of the Fock state scenario. 

On the other hand, if multiple photons are entering 
the system at different inputs, in general it is not possi- 
ble to simulate the system using coherent light. A sim- 
ple counter-example illustrates this. Consider the most 
trivial two mode quantum walk - a single beamsplitter 
with two inputs and two outputs. Should we put a single 
photon into each beamsplitter port the system under- 
goes the transformation |11) (|20) + \02))/V2. Thus 
upon performing number resolving detection, we will al- 
ways measure both photons at one port or both photons 
at the other (each with 50% probability), whereas one 
photon at each output will never be measured. With two 
coherent state inputs, there is no choice of phase rela- 
tionships which can mimic these statistics. The deflning 
difference here is that with single photons non-classical 
(HOM) interference is taking place. 

More formally we can consider the general case where 
an incident coherent state \a) propagates through an TV- 
port linear network consisting only of beamsplitters and 
phase-shifters. Such a network transforms the incident 
state into a product state of different coherent states. 



(19) 



Upon measuring single photons at output modes m and 
n, and no photons at all other modes, we obtain 



IV^OUt) 



A^£n4ne""''''^'|0) (20) 

/(m) fin) ^ ^ y 



We have presented a framework for describing multi- 
walker, discrete time quantum walks on arbitrary graphs, 
considering the cases of both distinguishable and indistin- 
guishable walkers. We presented several example appli- 
cations for this framework on different graph structures 
and discussed the speciflc example of photonic quantum 
walks. We discussed measurement and characterization 
of quantum walks and the issue of measurement corre- 
lations, which differentiate distinguishable versus indis- 
tinguishable walkers. We have shown that single walker 
walks have linear scaling in terms of the number of pa- 
rameters describing the system and that the complexity 
of the system scales exponentially with the number of 
walkers. Quantum walks have been shown to be univer- 
sal for quantum computation, and the framework we have 
presented provides the generic tools necessary to describe 
this. 
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